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ABSTRACT 

Multi-wavelength observations from radio to soft X-ray bands of large-scale 
galactic spiral structures offer synthesized and comprehensive views of nearby disk 
galaxies. In the presence of a massive dark-matter halo, the density-wave dynam- 
ics on galactic scales involves the stellar disk, the gas disk of interstellar medium 
(ISM), the magnetic field, and the cosmic-ray gas (CRG). In this paper, we ex- 
plore the dynamic and electromagnetic interplay between the magnetized ISM 
disk and CRG disk so that structural and diagnostic features of optical, infrared, 
and synchrotron radio-continuum emissions from a spiral galaxy can be physically 
understood. On timescales of galactic density waves, cosmic rays collectively may 
be treated as a relativistically hot tenuous gas fluid that is tied to the large-scale 
mean magnetic field in transverse bulk motions but moves otherwise differently 
along the magnetic field relative to the ISM. For both fast and slow magneto- 
hydrodynamic (MHD) density waves in a composite disk system of magnetized 
ISM and CRG, the minute CRG mass density enhancement is phase shifted rela- 
tive to the enhancement of parallel magnetic field. Owing to an extremely small 
number of cosmic rays, the large-scale magnetic field enhancement dominates in 
synchrotron radio-continuum emissions (as if the CRG is almost unperturbed) for 
spiral structural manifestations. In addition to the fast and slow MHD density 
waves, there also exists a suprathermal MHD wave mode by which CRG adjusts 
itself with an effective suprathermal sound speed close to the speed of light c. 

Key words: galaxies: cosmic rays — density waves — radio polarization — grav- 
itation — magnetic fields — interstellar medium 
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1. INTRODUCTION 

The entire system of a spiral galaxy consists of a massive dark-matter halo, an 
older stellar halo, a luminous stellar disk with an central ellipsoidal bulge, a thin 
disk of interstellar medium (ISM) consisting of a thermal gas plus dusts, relativis- 
tic cosmic rays, and a large-scale magnetic field (e.g., Woltjer 1965). The rotation 
curve and speed of a disk galaxy are dominantly controlled by the total mass 
(dark matter included) distribution in the system via gravity (e.g., Kormendy & 
Norman 1979; Kent 1986, 1987, 1988). The spiral structure in the disk seen in red 
light from the relatively old stellar population is usually broad and smooth (e.g. 
Elmegreen 1981), while the spiral arm structure in the disk seen in blue light from 
the young stellar population is typically brilliant, narrow, and sharp. Often, there 
exist dark narrow dust lanes lying along inner edges of the luminous spiral arms 
outlined by young O, B stars and Hll complexes. The large-scale spiral structure 
seen in total synchrotron radio-continuum emissions well tracks the optical spiral 
structure (e.g., the "Whirlpool galaxy" M51 [NGC 5194]; Mathewson, van der 
Kruit, & Brouw 1972; Neininger 1992; Berkhuijsen et al. 1997). Polarized radio- 
continuum emission arms also follow optical spiral arms as seen in several nearby 
spiral galaxies (e.g. M51, M31, NGC 2997). However, the spiral structure seen in 
polarized synchrotron radio-continuum emissions can, in the case of the nearby 
spiral galaxy NGC 6946 (Beck & Hoernes 1996; Fan & Lou 1996, 1999; Lou & Fan 
1998a, 2002; Ferguson et al. 1998; Frick et al. 2000; Lou 2002), appear interlaced 
with the optical spiral structure.^ Synchrotron radio-continuum emissions are pro- 
duced by relativistic cosmic-ray electrons gyrating around magnetic fields. The 
large-scale interrelations between optical and radio-continuum spiral structures 
seen in spiral galaxies hint at an intricate magnctoliydrodynamic (MHD) cou- 
pling among various seemingly unrelated processes on distinctly different scales 
(Lou & Fan 1997, 1998b, 2000a, b; Lou 2002). 

In addition to giving rise to a more or less flat disk rotation curve, the massive 
dark-matter halo also helps prevent rapid development of bar-type instabilities 
that have been shown to exist through numerical and theoretical studies (Miller, 
Prendergast, & Quirk 1970; Hohl 1971; Ostriker & Peebles 1973; Bardeen 1975; 
Shu et al. 2000). By joint effects of differential rotation, epicyclic oscillations, and 
self-gravity, such a disk system is still vulnerable to various instabilities (Safronov 
1960; Toomre 1964; Goldreich & Lynden-Bell 1965; Julian & Toomre 1966; Jog 
& Solomon 1984a, b; Binney & Tremaine 1987; Bertin et al. 1989a, b; Bertin 
& Lin 1996; Montenegro et al. 1999; Lou, Yuan, & Fan 2001). In the case 
of a stellar disk, such instabilities may ultimately lead to an effective increase of 
stellar velocity dispersion. In the case of a rotating gas disk, such instabilities may 
trigger cloud and star formations on global scales (Quirk 1972; Jog & Solomon 
1984a, b; Kennicutt 1989; Kennicutt ct al. 1994; Wang & Silk 1994; Elmegreen 
1995; Jog 1992, 1996; Silk 1997; Lou & Fan 1998b). As clouds and stars directly 

^ Two other late-type gas-rich spiral galaxies nearby, IC342 (Krause 1993; Crosth- 
watte et al. 2000) and M83 or NGC 5236 (Sukumar & Allen 1989; Neininger, Beck, 
Sukumar, & Allen 1993), also bear somewhat similar interlaced arm features of 
NGC 6946. 



4 



form in the gas disk of interstellar medium (ISM) and the magnetic and thermal 
energy densities in the ISM are roughly comparable, magnetic field should play an 
important role in afi'ecting the global star formation rate. In particular, magnetic 
field introduces additional MHD instabilities (Lou & Fan 1997, 1998a, 2000a, 
b; Lou et al. 2001a) that are independent of Toomre-type ring instabilities for 
axisymmetric disturbances and that can lead to nonaxisymmetric gravitational 
collapses (Lou 1996a; Lou et al. 2001a). 

Over the past decade, the advent of Infrared Astronomical Satellite (IRAS) 
and Infrared Space Observatory (ISO) has opened up the infrared window for sys- 
tematic, extensive, and unprecedented galactic observations. One important dis- 
covery of IRAS was the remarkably tight global correlation between integrated far- 
infrared and radio-continuum emissions from spiral galaxies (Dickey & Salpeter 
1984; Helou, Soifer Sz Rowan-Robinson 1985; de Jong et al. 1985; Wunderlich, 
Klein & Wielebinski 1987; Bicay & Helou 1990; Hclou 1991). This empirical fact 
implies a close physical connection between two apparently unrelated physical 
mechanisms, namely, thermal emissions from UV photon heated dusts in the ISM 
and synchrotron radio-continuum emissions from relativistic cosmic-ray electrons 
gyrating around galactic magnetic field. It was proposed that massive star for- 
mation in the ISM is likely the key process that somehow links these disparate 
radiative phenomena (Harwit & Pacini 1975; Helou et al. 1985; Condon, Ander- 
son, & Helou 1991). Recent high-resolution ISO observations of the nearby spiral 
galaxes NGC 6946 (e.g. Tuffs et al. 1996; Lu et al. 1996; Helou et al. 1996; 
Malhotra et al. 1996) and M31 (see extensive references in Berkhuijsen, Beck & 
Walterbos 2000) further revealed detailed spatial correlations and correspondences 
between large-scale infrared and radio structures as already partially evidenced 
in an earlier IRAS survey of spiral galaxies (Bicay & Helou 1990). We emphasize 
the importance of multi-wavelength observations of spiral galaxies (Lou, Walsh, 
Han & Fan 2002), because we believe that large-scale spiral MHD density waves 
are the underlying dynamic process that links and organizes the various ISM com- 
ponents through magnetic field and that gives rise to large-scale interrelations of 
spiral structures revealed in various wavebands (Lou & Fan 2000a, b). 

The familiar hydrodynamic density wave theory (Lin & Shu 1964, 1966) dealt 
primarily with large-scale perturbations in a thin stellar disk, using either a for- 
malism of stellar distribution function (e.g., Toomre 1964; Julian & Toomre 1966) 
or a fluid approach (e.g., Lin 1967b, 1987; Toomre 1977; Binney & Tremaine 1987; 
Bertin & Lin 1996). The inclusion of a less massive ISM disk was, in the past, 
regarded as a passive effect. Nevertheless, there has been growing evidence for the 
important dynamic as well as diagnostic roles of the ISM (e.g.. Jog & Solomon 
1984a, b). While smaller than the kinetic energy density of the galactic disk rota- 
tion, the energy densities of thermal ISM, magnetic field, and cosmic rays are com- 
parable (e.g., Lin 1967a). The increase of dynamical freedoms in the magnetized 
ISM may lead to different large-scale structural interrelations in multi-wavelength 
observations and the presence of magnetic field would significantly influence the 
global star formation rate in a nontrivial manner. Therefore, in order to under- 
stand large-scale infrared and synchrotron radio-continuum structures of spiral 
galaxies as well as their interrelations with the optical spiral structures, it be- 
comes necessary to take into account the effects of thermal ISM, magnetic field, 
and cosmic rays within the overall MHD density wave scenario. While the gravita- 
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tional effect dominates tlie dynamics of the massive stellar disk and of the mutual 
interaction between the stellar and ISM disks (Lou & Fan 1997, 1998b, 2000a, 
b), physical consequences of the dynamical interaction among the sub-systems 
remain to be specifically and thoroughly examined. 

While the ultimate goal of building up such an edifice remains challenging, it 
is nonetheless physically informative to extract necessary ingredients separately 
from such a comprehensive scenario in order to work out simpler yet non-trivial 
physical partial problems in concrete terms. This research strategy requires a 
systematic investigation on a series of partial problems involving various simplifi- 
cations. Meanwhile, one must clearly bear in mind the utility as well as limitations 
of the theoretical results thus derived. We hope to gain physical insights for the 
overall problem such that various aspects can be modeled quantitatively to con- 
front multi- wavelength observations of spiral galaxies (Lou et al. 2002). 

We have shown in earlier publications the possible existence of fast and slow 
MHD density waves^ in a differentially rotating thin gas disk embedded with an 
azimuthal magnetic field (Fan & Lou 1996; Lou & Fan 1998a). In particular, 
perturbation enhancements of thermal gas density and parallel magnetic field are 
roughly in phase for fast MHD density waves but are significantly phase shifted 
for slow MHD density waves with a phase difference > 7r/2. The features of fast 
MHD density waves may explain why optical and synchrotron radio-continuum 
structures nearly coincide in several nearby spiral galaxies (e.g. M51, M31, NGC 
2997 etc.), whereas the features of slow MHD density waves may explain why 
magnetic spiral arms lie in between the optical spiral arms in the galaxy NGC 
6946 and why such an interlaced spiral structure persists over a large spatial scale 
(Beck & Hoernes 1996; Fan & Lou 1996; Ferguson et al. 1998; Frick et al. 2000; 
Lou & Fan 1998a, 2002). 

An important new realm of applications for spiral MHD density wave theory 
has been developed recently for circumnuclear spiral arms and starburst "rings" 
on kiloparsec scales in the central regions around nuclei of disk galaxies such as 
NGC 1097, NGC 6951, NGC 2997, and NGC 2207 etc. (Lou, Yuan, Fan, & Leon 
2001). More specifically, when coupled with proper wave damping mechanisms in 
a circumnuclear disk, spiral MHD density waves play the key role of persistently 
removing angular momentum from the central magnetized gas disk in the cir- 
cumnuclear region, and the incessant accumulation of gas materials and magnetic 
flux in a roughly circular zone is eventually vulnerable to gravitational instabil- 
ities (e.g., Elmegreen 1994; Lou 1996a) and thus leads to the appearance of a 
starburst "ring" somewhere inside the modified inner Lindblad resonance (Lou 
et al. 2001b). At the current stage of observations, wavelet analyses and recon- 
structions on high-resolution optical images from Hubble Space Telescope (HST) 
have revealed circumnuclear spiral structures with or without nuclear bars in sev- 
eral dozens of spiral galaxies, including NGC 2207 (Yuan, private communications 
2001). The starburst galaxy M83 (NGC 5236) has been revealed to possess double 
circumnuclear ring and minibar on scales of several hundred parsecs (Elmegreen 



^ Note that incompressible Alfvenic fluctuations perpendicular to a magnetized 
rotating gas disk may exist (sort of bending MHD waves), but they do not directly 
couple to gravitational potential perturbations in the absence of mass density 
fluctuations (see the formulation contained in §2 and in Appendix C). 
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et al. 1998). In comparison, preliminary information of circumnuclear magnetic 
fields are only known recently for NGC 1097 (Beck et al. 1999) and barely for 
NGC 2997 (Han et al. 1999; Lou et al. 1999). Here again, it is the relativistic 
cosmic-ray electrons abundant in the circumnuclear environs that reveal valuable 
clues of trailing swirl patterns of circumnuclear magnetic fields. 

The central theme of this paper is to investigate the coupling of spiral MED 
density waves in the magnetized thermal ISM disk and the cosmic-ray gas (CRG). 
We are mainly motivated by the fact that the energy densities of thermal gas, 
magnetic field, and CRG are comparable (all on the order of ~ 10"^^ erg cm~^), 
the fact that galactic synchrotron radio-continuum emissions result from interac- 
tions of relativistic cosmic-ray electrons with the magnetic field, and the fact that 
magnetic field appears to be the underlying cause for the large-scale structural 
correlations as revealed by multi-wavelength observations of nearby spiral galax- 
ies. We do not include the effect of the massive stellar disk in the formulation 
merely for the sake of simplicity; yet we have enough confidence, based on our 
earlier analyses (Lou & Fan 1997, 1998b, 2000a, b), that the gross properties of 
fast and slow MHD density waves should remain in the magnetized thermal gas 
disk of ISM even when the more massive stellar disk is included (see Lou et al. 
2001a; Lou & Shen 2003). 

Our basic theoretical formalism of the interaction between the magnetized 
thermal gas and CRG closely follows that pioneered by Parker (1965, 1967, 1969) 
nearly four decades ago. The key starting point is to treat galactic cosmic rays 
collectively as a tenuous relativistically hot gas fluid for dynamic phenomena of 
large spatial scales (i.e. ^ gyroradii of cosmic rays) and low frequencies. In 
this paper, we therefore ignore resistive effects such that magnetic field is frozen 
into both ISM and CRG. The physical meaning of frozen-in condition is that 
bulk materials cannot go across magnetic field lines but are allowed to move 
along magnetic field lines. As ISM and CRG are two qualitatively different types 
of conducting fiuids, their bulk motions along magnetic field lines are different 
following the proper dynamics. In other words, both thermal ISM and CRG are 
tied to the galactic magnetic field in transverse bulk motions but are allowed to 
move relative to each other along magnetic field lines. By this analysis, it is then 
possible to assess the relative importance of the enhancements of parallel magnetic 
field and CRG mass density in synchrotron radio-continuum structures of spiral 
galaxies. 

The plan of the paper is as follows. Physical considerations and perturbation 
formulation of the problem are developed in §2. Dispersion relations of fast and 
slow MHD density waves and of the suprathermal MHD mode are derived in 
§3 using the tight-winding approximation. Galactic applications are described 
in §4. Notes and discussions are contained in §5. Appendix A contains some 
mathematical details of the key perturbation equation (3.12) for the convenience 
of reference and discussion. In Appendix B, we derive and summarize specific 
phase relationships of perturbation variables for the problem studied by Parker 
(1965) yet with the displacement current effect included. In Appendix C, we 
provide a first-principle formulation as the very basis of our MHD perturbation 
analysis. 
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2. FORMULATION OF THE PROBLEM 



For the mathematical description of the problem under consideration (see Ap- 
pendix C for more specifics), we adopt the cylindrical coordinate system {r,9,z) 
with the z- and disk rotation axes being coincident. For the stationary rotational 
equilibrium of the background, we presume the bulk azimuthal speed U0 of the 
CRG to be the same as the bulk azimuthal speed Vg of the thermal gas of ISM.^ 
The background magnetic field Bg is azimuthal to avoid the magnetic field winding 
dilemma due to the disk differential rotation (Lynden-Bell 1966; Roberts & Yuan 
1970). The profile of the background magneto-rotational equilibrium is character- 
ized by Ug = Vg = Qr and Bg = Fb/v, where 0(r) is the angular rotation rate and 
Fb is a constant. The stationary radial momentum equation for the background 
becomes 



dr 



f Bg a^TBg) CXPr 



(Lou 1996b; see Appendix C), where /Xo is the surface mass density of the thermal 
gas. Pa is the two-dimensional thermal gas pressure, (f)T is the total negative grav- 
itational potential, = So -|-rPo/[(r — l)c^], So is the vertically integrated mass 
density of the CRG, Po is the two-dimensional effective CRG pressure taken to be 
isotropic (Parker 1969), c is the speed of light, and P is the polytropic index of 
the CRG (e.g., P = 4/3 in a relativistically hot gas). The axisymmetric Poisson 
equation is 

with (j)T felt at the magnetized gas disk being determined by the distribution of 
the total mass density in the entire system. Within the solar system, the CRG 
was observed to be remarkably isotropic with an upper bound on anisotropy of 
~ 10"^ (Greisen 1960). In general, small-scale plasma instabilities are expected 
to reduce the anisotropy of the CRG to be < 1% (Lerche & Parker 1966; Lerche 
1967). Thermal gas pressure, effective CRG pressure, and magnetic energy density 
Bg/{8n) are roughly comparable and are on the order of ~ 10~^^ erg cm^'\ 

The thickness of the composite disk system is expediently taken to be infinitely 
thin here. In reality, the vertical quasi-static balance involves the confinement of 
the thermal ISM, CRG, and magnetic field under the action of gravity towards 
the galactic disk plane at 2: = from both sides. In the background rotational 
equilibirum, the radial Lorentz force in equation (2.1) vanishes (i.e., force-free) as 
rBe is taken to be a constant Fb. By this choice of Be that scales as r~\ one 
needs to invoke some processes in the vicinity of the galactic centre to avoid the 
outright singularity of Bg there. As the massive dark-matter halo contributes to 
Pt and thus to (^t, the disk rotation curve for the velocity Ves{r) of the stellar 
disk is actually used to infer (f)T from observations. Since the thermal energy 



^ In general, Vg and Ug may be different. In galactic contexts, this difference 
should be small, because the total gravity (that of dark matter halo included) 
dominates over pressure and Lorentz forces. Furthermore, Kelvin-Helmholtz in- 
stabilities induced by small differences in Vg and Ug may sustain a certain level of 
microturbulence that tends to smooth out such velocity shear. 
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density of the ISM is much smaller than the kinetic energy density of the galactic 
rotation and the magnetic field is taken to be force-free, the difference between 
the gas rotation velocity Ve(r) and the stellar rotation velocity Vgs{r) should be 
very small. In fact, one usually prescribes a V0{r) a priori based on observational 
input to construct an approximate background profile. 

Given a background profile for the thin rotating gas disk of magnetized ISM 
together with the CRG, it is fairly straightforward to write down MHD equa- 
tions for large-scale coplanar perturbations, that is, we limit the consideration to 
two-dimensional propagations of MHD wave perturbations tangential to the disk 
plane at z — 0. By this requirement, incompressible Alfvenic fluctuations involv- 
ing velocity and magnetic field perturbations perpendicular to the disk plane are 
excluded (Lou & Fan 1998a). Likewise, Parker instability (Parker 1966) and mag- 
netorotational instabilities (MRI; Balbus & Hawley 1998; Kim & Ostriker 2000) 
are also excluded. In the context of spiral MHD density waves, we are mainly 
interested in compressible fluctuations, for which the relevant MHD perturbation 
equations are given below explicitly. The radial component of the magnetic in- 
duction equation for the radial magnetic field perturbation is 

db, _ 1 djVeb,) 1 d{BeVr) . . 

dt r 89 r 09 ' ^ ' ^ 

where is the bulk radial thermal gas velocity perturbation, and the azimuthal 
magnetic induction equation for the azimuthal magnetic field perturbation bo is 

db^ ^ djVebr) _ d{BeVr) ,^ 
dt dr dr ' ^ ' ^ 

The divergence-free condition of the magnetic field perturbation b = {br, bg, 0) is 

i^ + l§ = 0. (2.5) 

r or r oti 

While the bulk radial velocity u,. of the CRG is the same as Vr as restricted by 
the frozen- in magnetic field on large scales, the bulk azimuthal fiow speed Ug of 
the CRG and the bulk azimuthal fiow speed Ve of the thermal gas can however be 
different along magnetic field lines. 

The radial component of the perturbed momentum equation is 



dVr _ 2lIoVgVg _ 2eoVgUg _^ 1^5^^ _ (// + €■) d{j)o -|- Pp) 

dt r(//o + eo) r{^io + eo) r 89 (//o + Co)^ dr 



1 d{p + P) f dzB, 



ifJ'o + Co) dr J 47rr(//o + Co) 



/ 



d{rbg) dbr 



dr 39 



dr 



(2.6) 



(Lou 1996b; see Appendix C), where e, p, and P are perturbations in //o, Co, 
Po, and Po, respectively, and is the negative self-gravity potential perturbation 
associated with the thermal ISM disk. The azimuthal momentum equation for 
the thermal ISM is 

dv0 ^ Vedve ^ l d{rVg) ^ ^ ^_dp ^ ld(f)_ 

dt r 89 r dr Ho'r 86 r 89 ' 
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and in parallel, the perturbed azimuthal momentum equation for the CRG is 

du^ ^ V^due_ ^ l d{rVe) ^ ^ _J_^ + 1^ (2 8) 

dt r 89 r dr e^r 86 r 86 

The mass conservation for the surface mass density perturbation // of the thermal 
ISM is 

8ii ^ 1 8{ii,rv,) ^ 8ve ^ Vedli _^ 
8t r 8r r 86 r 86 
The mass conservation for the mass density perturbation S of the CRG is 

as 1 dj^.ru^) S„ due , Vq dE _ 

dt^r dr ^ r 86^ r dd~^- ^^-^^^ 

The three-dimensional Poisson equation for the perturbed negative self-gravity 
potential (f) induced by fx becomes 

where the tiny contribution from the CRG to the negative self-gravity potential 
perturbation has been dropped, 6{z) is the Dirac delta function with argument z, 
and /Xo is related to the mass density of the thermal ISM by Po{r, z) = iio{r)d{z). 
For the magnetic field perturbation b, one of equations (2.3) and (2.4) may be 
spared when equation (2.5) is used instead. The perturbed polytropic relation for 
the two-dimensional ISM pressure perturbation p and the surface mass density 
perturbation p of the ISM becomes 

p = Clfi, (2.12) 

and similarly, the perturbed two-dimensional polytropic relation for the CRG is 

P = C^S , (2.13) 

where Cg and Cc are the polytropic sound speeds in the thermal ISM and in 
the CRG, respectively. In short, coplanar MHD density wave perturbations in 
the ISM disk and the bulk of the CRG are coupled electromagnetically via the 
large-scale mean magnetic field.* 

With the exp(za;t — im6) dependence implied for all MHD perturbation vari- 
ables in equations (2.3) — (2.13), where u) is the angular frequency in an inertial 
frame of reference and positive integer m indicates the number of spiral arms, one 
can reduce equations (2.3) — (2.13) to 

r{u> — mO)6r = —mBgVr , (2-14) 
'^^e = ^ , (2.15) 



^ Note that the large-scale dynamic coupling between the fluid stellar disk and 
the magnetized ISM gas disk is primarily gravitational (Jog & Solomon 1984a, b; 
Lou & Fan 1997, 1998b; Lou et al. 2001a). 
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d{rb,) 



m 



dr 



i{u3 — mVt)vr — 





dr (//o + €0) dr 



+ 



1 d{p + P) 



dzBa 



A'Kr{iJio + Co 



d{rhe) 
dr 



+ 



dr 



i{u} — mO,)ve + 



i{u! — mQ)ue + 



1 djrVo) 
r dr 

1 djrVe) 
r dr 



Vr- = — - 



■Ur — — - 



1 d{rfi,Vr) 
mil)iJL H ^— - 

r or 



%\LO — miZjS H 

r or 



imiic 



p_ 
P 

Co 

Vg = , 
U0 = O 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



The three-dimensional Poisson equation (2.11) leads to an integral representation 
of in terms of // (Shu 1970a; Lin & Lau 1979; Shu et al. 2000; GaUi et al. 2001; 
Lou & Fan 2001), which can be solved exactly for potential-density pairs in special 
cases (e.g., logarithmic spirals; Kalnajs 1971). In the present context, we assume 
a sufficiently short radial wavelength and thus invoke the tight-winding or WKBJ 
approximation. In this regime, it is possible to establish a local differential relation 
between (p and fi for the magnetized thermal gas disk of ISM (Lin & Shu 1964, 
1966; Shu 1970b). For spiral galaxies such as M51 and NGC 6946 of interest 
here, the tight- winding approximation is justifiable or, at least, should provide 
physically sensible results. We shall adopt this tight- winding approximation in 
our analytical analysis, as the solutions thus derived contain the information of 
MHD density waves in the presence of CRG fluid and of suprathermal MHD waves 
(see Parker 1965). 



3. DISPERSION RELATIONS OF SPIRAL MHD DENSITY WAVES 

Parker first investigated suprathermal hydromagnetic waves using the MHD 
perturbation equations in a thermal gas and some very hot suprathermal gas, such 
as interstellar cosmic rays, in a uniform magnetic field. The results of his pioneer 
analysis provide useful insights on our problem at hand, although the formula- 
tion here involves additional effects of disk differential rotation, self-gravity, and 
a curved magnetic field in the cylindrical geometry. Parker showed that, besides 
an incompressible transverse Alfven wave mode (see eq. [9] of Parker 1965), the 
conventional slow and fast hydromagnetic waves in the thermal gas are largely 
unaffected, except that there is a "hole" in the phase diagram for the fast mode 
in the direction perpendicular to the magnetic field (see Fig. 1 of Parker 1965; 
also see Appendix B here). The presence of this "hole" in the fast mode is related 
to the fact that bulk motions of thermal and suprathermal gases are tied to- 
gether electromagnetically in the direction transverse to the magnetic field. More 
importantly, there is an additional suprathermal mode, representing waves in the 
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suprathermal gas traveling effectively with the speed of sound in the superthermal 
gas alone (Appendix B). 

For spiral MHD density waves in a thin self-gravitating composite gas disk in 
rotation, we consider compressible MHD perturbations because mass density fluc- 
tuations are coupled to gravitational potential perturbation by Poisson's equation 
(2.11), so that a spiral gravitational potential field exerts a dynamical influence 
on compressible MHD perturbations. With a cylindrical geometry and in the 
presence of a thin differentially rotating disk, the analysis of MHD density waves 
becomes more tedious (Fan & Lou 1996; Lou & Fan 1998; Lou et al. 2001a), yet 
the analogs of fast and slow MHD waves involving a suprathermal gas as described 
by Parker (1965) should exist on the ground of physics. We expect that for fast 
MHD density waves, thermal gas density and magnetic field perturbations re- 
main more or less in phase in the tight-winding approximation, while thermal gas 
density and magnetic field perturbation enhancements remain significantly phase 
shifted for slow MHD density waves (a phase difference ^ 7r/2 but not exactly tt 
owing to the disk rotation). The existence of a somewhat modified suprathermal 
mode is thus anticipated as well. 

To derive the dispersion relations from equations (2.11) — (2.21) for the fast 
and slow spiral MHD density waves, we further take 



= $(r) exp 



i j k{s)df. 



(3.1) 



where k{r) is the radial wavenumber and $(r) is the slowly varying amplitude of 
the negative gravitational potential perturbation (j) (Goldreich & Tremaine 1978, 
1979). One then has 



dr 



ik{r)(j) -\- 



$'(r) 
$(r) 



(3.2) 



where the prime over $(r) denotes a radial derivative. In the tight-winding regime 
of a large /c, Poisson's equation (2.11) may be approximated by 



1 d{r 



1/2, 



^ 2nGi sgn(/c)/x + O 



z=0 



{kry 



(3.3) 



which contains a fractional error on the order of 0{kr)^^ (Shu 1970b; Goldreich 
& Tremaine 1979), where sgn(/c) = -|-1 for A; > and sgn(fc) = — 1 for k < 0. 
The cases of k > and k < correspond to leading and trailing spiral arms, 
respectively, and the integer m > gives the number of spiral arms. Equations 
(3.2) and (3.3) together give 



= 



AnriG sgn{k)fi 



2ikr + 2r$'/$ -|- 1 
A substitution of equation (3.4) back into equation (3.2) gives 

2r$'1 27riG' sgn{k)iJ, 



(3.4) 



dr 



2ikr + 



2ikr + 2r$'/$ + 1 



^1 1^ , 



(3.5) 
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which defines the complex coeflBcient T\. From equations (2.12), (2.18), (2.21) 
and (3.4), we can derive an equation of in terms of and dvr/dr 

irriK'^ HoVr/{2Qr) + i{uj — mQ)r~^ d{r j^i^v ,) / d.r 



^ {ijj- mny - m^Cyr^ + {m?/r^)4:-nriG sgn{k) fxj {2ikr + 2r$'/* + 1) 

= Av^ + B—^ , 
dr 

(3.6) 

which defines the two complex coefficients A and B. Prom equations (2.13) and 
(2.19) for perturbations in the CRG, one derives 

d{rVe) 



+ 



(3.7) 



r{uj — mQ) dr r{uj — mfl)eo 
From mass conservation (2.21) for perturbations of the CRG, one obtains 

?-rnE,/-rH,,/(2Q/-) + - inQ:)r-'d{vy.,u,) / dr 

~ {u - mfiy - rn'ClT.J{rH,) ' 

where n"^ = {2Q/r)d{rVg)/dr defines the epicyclic frequency k of the disk. From 
the mass conservation (2.20) for the thermal ISM, one has 

r{u3 — mVt)ii i d{rHoVr) 



(3.8) 



(3.9) 
(3.10) 

by using equation (3.10). A straightforward substitution of these equations into 
the radial momentum equation (2.17) leads to the following lengthy second-order 
ordinary differential equation (ODE) 



m//o m/ia dr 

The radial induction equation (2.14) may be arranged into the form of 

^ _ mBeVr 

r{u} — mVL) ' 

while the divergence-free condition (2.16) of h may be written as 



dr 



i{u) — mVt)vr 
20eo 



21]//o 



{l^o + Co) 



r{uj — mO) 



Av, + B 



dv. 



i d{riioVr) 



{jio + Co) \ r{uj - mO)eo 



TTi/Uo V dr J m//o dr 

imK^Uj./ (20r) + i{uj — mVt)d{rT^aUr) / dr / (rEo) 



[u - mVty - m^Q E„/(r2e,) 



+ 



20(u; - mfi) 



irn'^C\Vr 



TAAv, + B 



dVr 

dr 



r^{uj — mil) 

1 d 
(/io + eo)dr 



d ( d 
— < r — 
dr I dr 



CliAvr + B 



r{uj — mO) 

dVr 

dr 



1 d \ Cg,E4//n/-.-- »,/(2n;0 + - /n.n)rf(;-Eo»,) /(//■/(/■Ep)] 
(//o + Co) dr \ 



+ AVr + B 



+ 



d,Vr 



1 



(a; - m^Y - m^ClT.,/ {r^e,) 
d{po + Po) 



dr J (/io -I- dr 
rC5Eo/[(r - l)c^] d{po + Po) [imK^Ur/{2nr) + i{uj - mVl)d{rT.,Ur) / dr / {rT.,)] 



il^o + eoY 



dr 



{uj - mQY - m2C2Eo/(r2eo) 



(3.12) 
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where = / dzBl /[An^iio + to)] defines the Alfven speed Ca in the composite 
disk, and the complex coefiicients JF^, A, and B have aheady been defined earher 
by equations (3.5) and (3.6) (see Appendix A for details). With Vr = Ur, equation 
(3.12) contains the information of fast and slow spiral MHD density waves as well 
as the suprathermal MHD mode. 



3.1 SLOW SPIRAL MHD DENSITY WAVES 



We have previously (Fan & Lou 1996; Lou & Fan 1998a) derived, in the tight- 
winding approximation, the dispersion relation for slow MHD density waves under 
the conditions that Q = Vg/r ~ constant, uj — mQ. ~ mCA/T, Cs ~ C^, ^/Ca ~ /c, 
and kr ^ 1. Here, we require in addition that Cc 3> Cs and So//Uo -C 1. On the 
basis of Parker's analysis (Parker 1965), we expect that the basic features of slow 
MHD density waves should remain more or less intact and the associated CRG 
density fluctuation should be fairly weak. We set out to confirm these expectations 
by an analysis of equation (3.12) (see Appendix A) using the assumptions stated 
above. 

To the leading order of the tight- winding approximation with kr ^ 1 for slow 
MHD density waves, the three complex coefficients Ci, and C3 in equation {Al) 
are approximately given by 



+ 



i{u — mO) 



(a; — mfi) {jio + Co) (a; — m^iy — m^A/r^ 

i{oj - mn)C^So 



+ 



(/Xo + e.)[(a; - mny - m2QEo/(r2eo)] 



20;Uo 



(/Xo + Co) 



ir{uj — mO)' 



m[(u; — mfiy — m^A/r^ 
2zmf)C^So 



m 



— ik 



+ 



\k\ {u;-mny -m^A/r^ 



+ 



irriK^ / (20r) 



and 



C3 = 



(//o + Co) (a; — m^Y — m^A/r^ 
i{u) — mVt)^'^ 



+ 



[to 



uiVlY — m? A / r"^ (yUo + Co) 
irriK? / (20r) 



IK 



2VL{uj - mVL) 



r{u — mfl)€o (w — mfiy — m^C^So/ (r^eo) 
27rGiJo imKy{2nr) 



ik 



\k\ {u) — mVty — m^A/r'^ 



(3.13) 



(3.14) 



(3.15) 
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separately with A = C| — 27rG'//o/|/^| (see Appendix A). Ignoring higher-order 
terms, the dispersion relation of slow MHD density waves comes out to be 



. ^ (Ato + Co) (//o + Co) |A;| 



(3.16) 



.(/"o + eo) C^Eo(/io + eo 

Because //o S> or So, the slow MHD density wave is indeed not very much 
affected by the presence of a tenuous CRG (Fan & Lou 1996; Lou & Fan 1998a, 
2001). As expected, perturbation enhancements of surface mass density n of the 
thermal ISM gas and parallel magnetic field bg are significantly phase shifted with 
a phase difference > tt/2. The particular aspect of interest here is the phase 
relation and magnitude of CRG density fluctuation E relative to thermal gas 
density fluctuation fj,. By equations (3.6) and (3.8), we have for large kr 



fi/fio (a; — mQ)^ — m^C^Eo/(r^eo) 



(3.17) 



For neutral or stable slow MHD density waves with A > 0, the inequality [co — 
mfiy < m^A/r^ -C m^C^Eo/(r^eo) holds such that E/E^ and /////o are in phase 
with 

— ~ ^ . (3 18) 

Eq C^Eq /io 

As Co and Eo are on the same order of magnitude and A ^ C^, the relative 
fluctuation E/Eo of CRG density is much less than the relative fluctuation /x///o 
of thermal gas density. Therefore, the main contribution to the large-scale vari- 
ation of synchrotron radio-continuum emission comes from the enhancement of 
magnetic held. In other words, the signiflcant phase difference ( > tt/2) between 
the enhancements of CRG density and parallel magnetic fleld would not reduce 
the contrast of synchrotron radio-continuum emission associated with slow MHD 
density waves in a signiflcant manner. 

As the presence of tenuous CRG does not influence the propagation of slow 
MHD density waves in a signiflcant manner, it follows that the dispersion relation 
for slow MHD density waves in a disk with a strong differential rotation should 
remain more or less the same as derived earlier (Lou & Fan 1998a). 



3.2 FAST SPIRAL MHD DENSITY WAVES 



To derive the dispersion relation of fast MHD density waves in the presence 
of a CRG, we take that {to - mQ) - (C| + ClY^^k, Cg - Ca, ^/Ca - fc, 
kr ^ 1, Cc ^ Cs, Eo <^ //o, and m/{kr) » Cs/Cc-, where the last inequality 
guarantees that the direction of wave propagation be sufficiently away from the 
radial direction. It is then straightforward to show (see Appendix A) that 

^1 ~ 7 7^ 



i{u - mn)C^Eo ^ ^ 

^ (Ato + - mO)2 - m2C2Eo/(r2eo)] ' 
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where the third term on the right-hand side due to the CRG is small, 



\k\ {uj — mVlY — m? /S./r'^ 



(3.20) 



and 



C3 = i{u) 



mVL) — 

2n 



(a; — mf2)K^ 



(a*o + Co) (oJ - mi).y - 777,2 A/r^ _|_ ^^^^^ _ ^^^^ 

imK^/{2Vtr) + z(a; - mVt)d\n{rT.o) / dr 



(3.21) 



(//o + eo) r(a; - ml]) (a; - mVLY - m'ClT.,/{r''e,) 



in equation {Al). 
thus given by 



The dispersion relation of fast spiral MHD density waves is 



(a; — vn^iy 



iHo + Co 



Cgl-lo 



27rG//o 



(/^o + Co) \k\ 



(3.22) 



As expected, the presence of CRG does not affect very much the propagation of 
fast MHD density waves as long as m/{kr) 3> Cg/Cc- This inequality guarantees 
that a propagation of fast MHD density waves is sufficiently away from the radial 
direction that is perpendicular to the background magnetic field (see Appendix 
B). One can readily show that the perturbation enhancements of thermal gas 
density and parallel magnetic field are largely in phase to the leading order of 
large kr. From equation (3.17), it is clear that the perturbation enhancements of 
thermal ISM density and CRG density are out of phase (i.e., a phase difference 
of ~ tt). Furthermore, 

- ■ (3.23) 



which gives a much smaller S/Eo relative to /x//Xo. Nevertheless, in comparison 
to the case of slow MHD density waves, E/Eo is larger by a factor k^r^ jw? in 
the case of fast MHD density waves. Although perturbation enhancements of 
CRG density and parallel magnetic field are out of phase, the synchrotron radio- 
continuum emission is dominantly determined by the perturbation enhancement 
of large-scale parallel magnetic field. 



3.3 THE SUPRATHERMAL MODE 



To derive the dispersion relation of the suprathermal mode, we take that 
(a; - mQ) ~ mCc/r, Cg ~ Ca, ^/Ca k, Cc Cg, Eo < //o, and kr > 1. It 
follows that in equation (^1) 

^1 ~ 7 7^ 



i{u - mn)C^Eo ^ ' ^ 

^ (Ato + - mO)2 - m2C2Eo/(r2eo)] ' 
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k 



{u> — mO) 



and 



C3 = i{uj — mVL) 



\k\ {cu -mny -m^A/r^ ' 

i^o (a; — mO)^^ 



{Ho + 60) {lo — mQy — m^A/r^ 
(see Appendix A). By dropping several higher-order terms, we arrive at 



(3.25) 
(3.26) 



{uu — mfiy 



+ 



+ 



(/Uo + Co) 
2 

(/io + Co 



(/io + Co 



(A^o + Co) 



H — >{co - mil) 



(3.27) 



ilJ'o + Co) 



27rG/Uo 

|A;| 



9 — 'J ) 



which contains both the suprathermal wave mode and the fast MHD density 
wave mode. If the term {u — mn)^/c^C^So/(/Xo + Co) were absent, then relation 
(3.27) may be readily expressed as a multiplication of two simple factors. To the 
present level of approximation, the dispersion relation for the suprathermal mode 
is approximately given by 



+ 



(Ato + Co) |A;| 



(3.28) 



r^Co (/io + Co) (//o + Co) 
Again, from the relation (3.17), it is seen that E/So and /////o are in phase with 

(3.29) 



where the magnitude of E/So may be comparable to or larger than that of (JLo 
for the suprathermal mode. 

The velocity perturbations of the suprathermal mode can be derived from the 
following relationship 



Ue 



mC^Eo imn'^Ur/ (2Vtr) + i{uj — mQ)[urd\n{rTio) / dr + dur/dr] 



(w — mVty — m^C^So/ (r^eo) 



+ 



• 2 



(3.30) 



For an order-of- magnitude estimates, we have 



krCgSo 



(3.31) 



For ttr ~ 1 km s~\ the azimuthal velocity perturbation ue could be as large as 
~ 10'' km s~\ For Ug > C^, one should be seriously concerned with streaming 
instabilities along the magnetic field. Parker (1969) pointed out the possibility of 
Landau damping for the suprathermal mode in the CRG. 
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4. APPLICATIONS TO SPIRAL GALAXIES 

To apply the results of our analyses to large-scale spiral structures of gas- 
rich disk galaxies seen in synchrotron radio-continuum emissions, one needs to 
have an idea for the plausible state of ORG in a typical magnetized disk galaxy. 
If the background CRG as a whole (either a fat disk or an oblated spheroid) 
largely corotates with the stellar and magnetized gas disks while retains large- 
scale axisymmetry, then our theoretical results indicate the following. Firstly, in 
the presence CRG density fluctuations on large spatial scales, the suprathermal 
mode propagates essentially with the cosmic-ray sound speed (always less than but 
very close to the speed of light c) more or less along magnetic field lines to smooth 
out CRG fluctuations. In this case, the large-scale distribution of CRG is likely 
to vary in the radial direction due to the presence of a mean azimuthal galactic 
magnetic field. Secondly, CRG density fluctuations associated with slow MHD 
density waves are extremely weak as compared to thermal gas density fluctuations. 
While the perturbation enhancements of CRG density and parallel magnetic held 
are signiflcantly phase shifted (i.e., > 7r/2), the enhancement in synchrotron 
radio-continuum emission is mainly determined by the enhancement of parallel 
magnetic fleld. Thirdly, CRG density fluctuations associated with fast MHD 
density waves are also sufficiently weaker than thermal gas density fluctuations to 
warrant that the enhancement of synchrotron radio-continuum emission basically 
goes with the enhancement of parallel magnetic fleld; the latter anti-correlates 
with the minute enhancement of CRG density fluctuation. 

For an actual spiral galaxy, several complications and uncertainties arise. The 
key question ties to the actual distribution of relativistically hot CRG in a typical 
disk spiral galaxy. We further elaborate several relevant aspects below. 

1. It is known observationally that the total luminosities of infrared and 
nonthermal radio-continuum emissions from spiral galaxies are tightly correlated 
(Dickey & Salpeter 1984; Helou et al. 1985; Helou & Bicay 1993). This empirical 
fact leads to the important notion that somehow the formation of massive stars 
and the production of relativistic cosmic rays are correlated through a chain of 
complex physical processes yet to be explored and understood. For example, an 
active region of massive star formation may produce, perhaps in the statistical 
sense, more explosions of type II supernovae that are favored sources of cosmic ray 
production (up to energy of ~ 10^'' eV). It is also plausible that MHD termination 
shocks conflning fast stellar winds from numerous young O, B stars sustain the 
supply of galactic cosmic rays in a signiflcant manner. It thus seems sensible that 
on large spatial scales, optical spiral gas arms of relatively high densities should 
contain more sources for cosmic ray production. 

2. Interstellar cloud complexes, where young massive stars are continuously 
forming on smaller scales, are generally thought to be initially triggered by Parker 
instability (i.e., magnetic Ray leigh- Taylor instability, see Parker 1966, 1969; Shu 
1974; Mouschovias 1996) which involves downward slide of condensed thermal gas 
into magnetic valleys and upward expansions of magnetic flux tubes further in- 
flated by CRG in the dimension perpendicular to the galactic disk. To account for 
such processes, one must allow for MHD fluctuations perpendicular to a galactic 
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plane. The development of Parker instability disturbs galactic magnetic field on 
horizontal size scales of ~ 1 kpc or so. Subsequent activities of star formation can 
lead to intensified but disordered magnetic fields on smaller spatial scales via ISM 
turbulence. Therefore, the total magnetic field strength may increase with star 
formation activities in high-density gas arms. As enhanced magnetic fields in this 
manner tend to be disordered, the degree of polarization in radio-continuum emis- 
sions should be correspondingly low. Therefore for fast MHD density waves with 
in-phase gas density and magnetic field enhancements, manifestations of total and 
polarized radio-continuum emissions are in competition. That is, along a high- 
density gas arm with strong star formation activities, the total radio-continuum 
emission would be very strong whereas the degree of polarization would become 
weak. For moderate star formation activities along high-density gas arms, total 
emissions are generally stronger than polarized emissions but one may still detect 
degrees of polarizations. For both M51 (Neininger 1992) and NGC 2997 (Han et 
al. 1999), polarized radio intensity peaks closely follow the inner edges of optical 
spiral arms and coincide with narrow dark dust lanes. For M31 (Beck et al. 1980), 
magnetic "torus" coincide with the "ring" (for the controversy of "ring" or tight- 
winding spiral arms for M31, the reader is referred to references in Berkhuijsen et 
al. 2000). On the other hand, for slow MHD density waves with phase-shifted gas 
density and magnetic field enhancements, polarized radio intensity arms are much 
less affected by star formation activities along high-density gas arms where total 
radio emissions still remain to be strong. In terms of overall structural patterns, 
NGC 6946 fits this scenario weU (Lou & Fan 2002). 

3. Observations indicate that large-scale galactic magnetic field normally lies 
more or less in the disk plane of a spiral galaxy. Within the disk of a nearly face- 
on spiral galaxy (e.g., M51, NGC 2997 etc), polarized radio-continuum emissions 
reveal a global swirl-like pattern in the same sense of optical spiral arms (Neininger 
1992; Berkhuijsen et al. 1997; Han et al. 1999).^ The usual impression that 
magnetic field remains aligned and connected along spiral arms would be difficult 
to reconcile with the existence of a strong disk differential rotation (Lynden- 
Bell 1966; Roberts & Yuan 1970). Disk differential rotation would wrap up a 
connected magnetic field line along a spiral arm in a few turns of galactic rotation. 
Similar to the winding dilemma of material spiral arms (Goldreich & Lynden- 
Bell 1965; Lin & Shu 1966), this winding dilemma for the magnetic field may 
be resolved by considering large-scale azimuthal galactic magnetic field strongly 
distorted by galactic MHD density waves in a systematic manner; it is plausible 
that enhanced magnetic field locally orients roughly along a spiral arm but is not 
actually connected along the spiral arm. In this scenario, magnetic spiral arms 
appear as a result of the collective behavior of phase-organized, distorted rings of 
magnetic field and gas stream lines®. We note that M51 seen in linearly polarized 
optical emissions shows a more circular or circumferential pattern of magnetic 
field configuration (e.g., Scarrott, Ward-Thompson & Warren-Smith 1987). 

4. For the scenario of enhanced magnetic field aligned and connected along a 

^ At longer radio wavelengths, the effect of Faraday rotation gives rise to more 
distorted or irregular swirl-like patterns (e.g., Berkhuijsen et al. 1997). 

® This is very much like the kinematic description of density waves in a stellar 
disk by visualizing the collective behavior of phase-organized, distorted rings of 
mean stellar orbits as suggested by Kalnajs (1973). 
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high-density spiral gas arm, cosmic rays produced by star formation and type II 
supernovae would be somewhat confined to retain relatively high number density 
along the arm. While for the scenario of distorted rings of magnetic field associated 
with MHD density waves, cosmic rays produced in high gas density spiral arms 
can redistribute themselves more readily along distorted rings of magnetic field, 
maintaining a large-scale radial gradient across rings of magnetic field at the same 
time. For a mature late-type spiral galaxy, the level of cosmic rays has probably 
reached the maximum that a galaxy can contain and remains more or less steady 
(i.e., leakage and decays are globally balanced by production of cosmic rays). In 
the azimuthal direction, a relatively high concentration of cosmic rays can be 
smoothed out via propagation of the suprathermal mode (Parker 1965, 1969), by 
various diffusion processes (see Cesarsky 1980) and by scattering processes. 

5. It is the more common situation that the spiral pattern of the total (ran- 
dom and organized magnetic fields plus cosmic rays) synchrotron radio-continuum 
emission more or less coincides with the optical spiral pattern in nearby spiral 
galaxies. Typically, a spiral arm of total synchrotron radio-continuum emission 
is broader and may extend further outside optical arms. This demonstrates the 
ability of cosmic rays to diffuse across magnetic field and to infiltrate a much 
larger spatial volume than their original sources of production. Qualitatively, one 
tends to interpret the optical-radio spiral pattern correlation in terms of some- 
what high concentration of cosmic rays and enhanced random magnetic fields 
along a spiral arm of high thermal gas density where stars are actively forming. 
As large-scale organized magnetic field may not be completely disrupted along 
star-forming regions, coherent structures in polarized radio-continuum emissions 
can be still detected. 

6. In addition to the fact that total luminosities of infrared and synchrotron 
radio-continuum emission from spiral galaxies are tightly correlated, new evidence 
is forthcoming for the infrared- radio correlation in spiral structural patterns within 
disk galaxies (e.g., M31, see references in Berkhuijsen et al. 2000 and Lou & Fan 
2000a, b), which is quite similar to optical-radio correlation in spiral patterns. 
With increasing sensitivity and angular resolutions of infrared observations, we 
expect more examples of detailed spiral arm correlations. This trend of infrared- 
radio correlation in spiral patterns would be consistent with the general perception 
of spiral MHD density waves, even though the chain of physical processes on 
various scales that lead to such correlated radiative manifestations remain to be 
thoroughly understood. 

7. In identifying fast MHD density waves with large-scale optical and radio 
spiral patterns in M51 (Fan & Lou 1996; Lou & Fan 1998a), we follow a few 
clues. Firstly, both polarized intensity maps and degree-of-polarization maps for 
synchrotron radio-continuum emissions from M51 (Neininger 1992; Berkhuijsen 
et al. 1997) reveal a global spiral pattern well correlated with the optical spiral 
pattern. This is the strongest evidence that large-scale galactic magnetic fields are 
involved in the dynamics of density waves. Secondly, both polarized intensity and 
degree of polarization are usually strong along optical spiral arms where thermal 
gas concentration is relatively high. The coincident enhancement of thermal gas 
density and magnetic field is expected from the basic property of fast MHD density 
waves. This implies that the process of star formation does not completely disrupt 
large-scale galactic magnetic field, although one cannot tell immediately whether 
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CRG density is much higher in spiral arms than in interarm regions. Even if CRG 
density distributes rather evenly in the azimuthal direction, it is still expected to 
see enhanced polarized and total radio-continuum intensities. Finally, from the 
perspective of self excitation and maintenance of MHD density waves'', fast MHD 
density waves are preferentially swing amplified in a disk of strong diff'erential 
rotation (Fan & Lou 1997) which is the case for M51 as weU as for NGC 2997. 

8. In identifying slow MHD density waves with large-scale optical and radio 
spiral patterns in the spiral galaxy NGC 6946 (Beck & Hoernes 1996; Beck et al. 
1996; Fan & Lou 1996; Ferguson et al. 1998; Frick et al. 2000; Lou & Fan 1998a, 
2002), we would like to clarify a few points. Firstly, although somewhat fuzzy, 
large-scale radio spiral pattern in NGC 6946 again indicates the involvement of 
magnetic field in global density waves. The total synchrotron radio-continuum 
emission along optical spiral arms is strong (see Fig. 3 of Beck & Hoernes 1996). 
This can be explained, as usual, in terms of the enhancement of small-scale random 
magnetic fields and cosmic ray production along zones of relatively high thermal 
gas concentration where stars are continuously borne. Secondly, intensity and 
degree of polarized radio emission are sufficiently strong to reveal global magnetic 
arms interlaced with optical spiral arms^ (see Fig. 2 of Beck & Hoernes 1996 
and Frick et al. 2000). For slow MHD density waves (Fan & Lou 1996; Lou & 
Fan 1998a), the enhancements of parallel magnetic field and thermal gas density 
are significantly phase shifted (i.e., > 't/2), and the associated fiuctuation of 
CRG density is so weak as can be neglected. It is therefore natural to interpret 
the prominent magnetic arms in polarized radio-continuum intensity in terms of 
enhancements of large-scale organized magnetic field, associated with slow MHD 
density waves, submerged in an omnipresent CRG. In this scenario, magnetic 
arms are interlaced with the optical spiral arms as a result of significant phase 
shift between enhancements of parallel magnetic field and thermal gas density. 
The fact that the polarized radio intensity is extremely low along optical arms 
(or total radio intensity arms; see Fig. 3 of Beck & Hoernes 1996) may be due to 
a nearly complete disruption of an already weakened large-scale magnetic field or 
due to an almost flat distribution of CRG along the azimuthal direction^ or both. 
At least, the CRG density cannot be too much higher across optical arms than in 
interarm regions; otherwise, the presence of a weak organized large-scale magnetic 
field may lead to a polarized radio-continuum intensity peak within optical arms 
as well. 

9. We have developed very recently (Lou & Fan 2002; Lou 2002) models for 
stationary fast and slow MHD density waves with logarithmic spirals in a mag- 



^ The presence of a companion (NGC 5195) around M51 (NGC 5194) is prob- 
ably responsible only for the distortion in the outer spiral pattern as a result of 
tidal interaction (Toomre & Toomre 1972; Elmegreen et al. 1989); the mner spiral 
pattern may be self-sustained (Rix & Rieke 1993; Lin 1996). 

* Even in the interarm regions, total radio-continuum emission is several times 
stronger than the peak polarized radio-continuum intensity. However, an interarm 
bump in total radio emission can be attributed to polarized intensity peak. 

^ That is, CRG density enhancement across star-forming optical arms may be 
smoothed out sufficiently fast such that the distribution remains more or less 
uniform in the azimuthal direction along magnetic field lines. 



21 



netized singular isothermal disk (MSID) with a flat rotation curve. The Poisson 
equation can be solved exactly without invoking the tight-winding approximation 
(Kalnajs 1971; Shu et al. 2000) so that the pitch angle of a logarithmic spiral 
can be arbitrary. Our main motivation is to address the issue of sustaining an 
extended slow MHD density wave pattern within a disk of a flat rotation curve 
(Sofue 1996; Ferguson et al. 1998; Frick et al. 2000). Given what we have learned 
from this analysis, it is possible to further include a CRG into the formalism of 
Lou & Fan (2002). It is expected that minute enhancements of logarithmic spiral 
arms in the CRG would be phase-shifted relative to enhancements of logarithmic 
spiral arms of magnetic field for either fast and slow MHD density waves. An 
additional stationary MSID configuration with logarithmic arms is also expected, 
corresponding to the suprathermal mode (Parker 1965). 

5. NOTES AND DISCUSSIONS 

We have studied fast and slow MHD density waves in a rotating magnetized 
thermal gas disk in the presence of CRG that is treated as a tenuous, relativisti- 
cally hot gas fluid for large-scale and low-frequency dynamic phenomena. In this 
formalism, thermal gas and CRG are tied to the galactic magnetic field in trans- 
verse bulk motions but are allowed to move relative to each other along magnetic 
field lines. One way to better appreciate fast and slow MHD density waves (Fan 
& Lou 1996; Lou & Fan 1998a; Lou et al. 2001a; Lou 2002) is to first under- 
stand the basic properties of fast and slow MHD waves in a uniformly magnetized 
medium. Similarly, the way to fully appreciate fast and slow MHD density waves 
in the presence of CRG in the current context is to acquaint oneself with the 
basic properties of MHD wave modes in a uniformly magnetized thermal medium 
submerged in a uniform CRG pioneered by Parker (Parker 1965, 1967, 1969; see 
also Appendix B for more specifics). 

The analysis of Parker (1965) leads to following basic theoretical facts. First, 
the conventional slow MHD wave in the thermal gas is not very much affected by 
the presence of CRG, and the associated CRG mass density fluctuation, which 
anti-correlates with the parallel magnetic field perturbation (Appendix B), is tiny 
compared to the thermal gas density fiuctuation. Secondly, the conventional fast 
MHD wave in the thermal gas, propagating sufficiently away from the direction 
perpendicular to the background magnetic field, is not very much affected by the 
presence of CRG, and the associated CRG mass density fluctuations, which also 
anti-correlates with the parallel magnetic fleld perturbation, is also tiny compared 
to the thermal gas density fluctuation. In the propagation direction perpendic- 
ular to the magnetic fleld, the conventional fast MHD wave is suppressed (see 
Fig. 1 of Parker 1965). And flnally, there exists a suprathermal mode resulting 
primarily from the compression of a relativistically hot CRG with a characteristic 
sound speed close to the speed of light c. Enhancements of CRG density, paral- 
lel magnetic fleld, and thermal gas density are all in phase for the suprathermal 
mode. 

The speciflc aims of performing this analysis are to derive the phase relation- 
ships between the perturbation enhancements of parallel magnetic fleld and CRG 
mass density associated with fast and slow MHD density waves, to conflrm that 
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the phase relationships between the perturbation enhancements of parallel mag- 
netic field and thermal gas density associated with fast and slow MHD density- 
waves in the presence of CRG remain more or less the same as in the case where 
CRG is absent (Fan & Lou 1996; Lou & Fan 1998a), and to assess the net effect 
on synchrotron radio-continuum emissions as a result of fluctuations in parallel 
magnetic field strength and in CRG mass density. 

The basic conclusions of our analysis are the following. (1) In the presence of 
CRG, fast and slow MHD density waves are only slightly modified, and there exists 
an additional generalized suprathermal mode. (2) As expected, the perturbation 
enhancements of parallel magnetic field and the thermal gas are in phase for fast 
MHD density waves and are significantly phase shifted (i.e. > 7r/2) for slow MHD 
density waves. (3) The perturbation enhancements of the thermal gas and CRG 
mass density are in phase for slow MHD density waves and are signifcantly phase 
shifted (i.e., > 7r/2) for fast MHD density waves; in other words, the perturbation 
enhancements of parallel magnetic field and CRG mass density are significantly 
phase shifted (i.e., > 7r/2) for both fast and slow MHD density waves. (4) As 
CRG mass density fluctuation associated fast and slow MHD density waves is 
tiny, the radio-continuum structures seen in synchrotron emissions are primarily 
determined by magnetic field enhancements. 

One restriction of our analysis is to require that velocity and magnetic field 
perturbations lie within the gas disk plane. This constraint, which leads to a 
great deal of simplification in the analysis, does exclude the possible occurrence 
of Parker instability (i.e., magnetic Ray leigh- Taylor instability) in the vertical 
direction across the magnetized disk (Parker 1966, 1967; Shu 1974; Mouschovias 
1996 and extensive references therein). The Parker instability gives rise to mag- 
netic loops or bubbles inflated by cosmic rays on both sides of the disk and causes 
condensation or conglomeration of thermal gas clouds into magnetic valleys on 
spatial scales in the order of ~ 1 kpc (give-and-take a factor of 2) along the mean 
magnetic held. As a result of vertical undulation of magnetic fleld lines and redis- 
tribution of CRG, synchrotron intensity can vary on the scale of ~ 1 kpc. Since 
fast and slow MHD density waves vary on spatial scales > several kpc, the over- 
all enhancement of synchrotron emission resulting from transverse compression 
of large-scale magnetic fleld within the disk will not be fundamentally altered by 
relatively small scale magnetic bubbles pertruding out of the disk plane. In short, 
even though quantitative calculations of synchrotron contrast must take conse- 
quences of Parker instability into account, the enhancement of large-scale galactic 
magnetic fleld plays a dominant role for increasing synchrotron radio-continuum 
intensity. 

Another restriction of our analysis is to require that waves propagate within 
the disk plane, with vertical wavenumber = 0. As noted recently in Lou et al. 
(2001a), this assumption excludes a class of magnetorotational instabilities (MRIs) 
given the specifled disk and fleld geometries (e.g., Balbus & Hawley 1992; Terquem 
& Papaloizou 1996; Kim & Ostriker 2000). The well-known axisymmetric MRIs 
(Velikhov 1959; Chandrasekhar 1960, 1961; Balbus & Hawley 1991, 1998) involve a 
disk differential rotation, a weak magnetic fleld component parallel to the rotation 
axis, and vertical perturbation variations with kz 7^ 0. For nonaxisymmetric MRIs 
in a disk with a toroidal magnetic fleld, it is essential to realize the fact k^ ^ that 
couples the perturbation equations for velocity and magnetic fleld perturbations 
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perpendicular to the disk and the equations for MHD perturbations coplanar with 
the disk (Lou et al. 2001a). 
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APPENDIX A 

By requiring Ur = fr, the lengthy equation (3.12) for Vr contains the infor- 
mation of fast and slow MHD density waves and the suprathermal mode (Parker 
1965) and can be cast into the compact mathematical form of 

^ d'^Vr ^ dVr ^ , . s 

Ci—i+C,-^+C3Vr = 0, {Al) 



where the three complex coefficients Ci, C2, and C3 are defined by 
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(A4) 

The complex coefficients ^i, ^, and B are defined by equations (3.5) and (3.6) in 
the main text. 



APPENDIX B 



For the convenience of reference and discussion, we derive and summarize 
here several pertinent phase relationships among perturbation variables, with a 
generalization of the original analysis by Parker (1965) for hydromagnetic pertur- 
bations in a uniformly magnetized medium of thermal and suprathermal gases. 
To accommodate the possible parameter regime of = i?/[47r(po + ^o)]^^^ ^ c 
with c being the speed of light, we include the effect of displacement current per- 
turbation for the Lorentz force term in the momentum equation. Otherwise, the 
Cartesian coordinates {x,y,z), notations, and definitions all remain essentially 
the same as those adopted by Parker (1965). 

From the mass conservation and 2;— component of the momentum equation for 
the thermal gas, one has 



p{uj^ — k'^a^ cos^ 0) 
Lopok sin 9 



(51) 
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where Uy is the bulk velocity perturbation of the thermal gas perpendicular to 
the background magnetic field Bz^ u is the angular frequency, k is the total 
wavenumber, 6 is the angle in the y^;— plane between the wave vector k and the 
background magnetic field Bz^ a is the sound speed of the thermal gas, po is 
the uniform background thermal gas density, and p is the thermal gas density 
perturbation. In parallel, from the mass conservation and 2;— component of the 
momentum equation for the suprathermal gas, one has 

5{u'' - k^h^ cos^ 6) 

where Vy is the bulk velocity perturbation of the suprathermal gas perpendicular 
to the background magnetic field Bz, b is the sound speed in the suprathermal gas, 
So is the uniform background suprathermal gas density, and 6 is the suprathermal 
gas density perturbation. By the requirement that the bulk velocities of thermal 
gas and superathermal gas are the same transverse to the background magnetic 
field Bz, namely Uy = Vy, we derive from equations (Bl) and {B2) 

6 p{u;^ — k'^a'^ cos^ 9) 

" p,(uj^ - k^¥ cos^ 0) ' ^ ' 

where S/So and p/po are the relative mass density fluctuations in the suprathermal 
and thermal gases, respectively. 

By combining the 2;— component of the magnetic induction equation 

u}b^ = BksmdUy {B4) 

with equation (Bl) where 6^ is the 2— component of the magnetic field pertur- 
bation, one has the relative fluctuation in parallel magnetic field b^/B and the 
relative mass density fluctuation in thermal gas p/po related by 

b^ _ p{u^ - Pg^ cos^ 6') 

B ~ p,u^ ' ^ ^ 

The incompressible transverse Alfven mode involves — and b^ ^ 0, and 
solution (9) of Parker (1965) for its dispersion relation now takes the form of 

" = (1 + wyc^) ' 

where W = B/[4n{pa + So)Y^'^- The modified dimensionless dispersion relation for 
the fast, slow, and suprathermal modes now becomes 

_ M') - cos' 6) {U' - n' cos' 9) 

U' sin' ^[(1 + a)U'-{l + a/n')n' cos' 9] _ ^ (S7) 



(l + a/?i2)(l + M^7c2) 



= 0, 



where U = uj/{ka), M' = W'/[a'{l + W'/c')], n = b/a, a = 5,b'/{p,a') (for 
a comparison, see Parker's eqn [15] and the relevant definitions). The factor 
(1 + W'/c') that appears in various places is due to the displacement current 
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effect. We now examine systematically the phase relationships among various 
perturbation enhancements such as p/po, b^/B, and 5/5o for the fast, slow, and 
suprathermal modes separately, using the phase relations (S3) and (S5). 

1. The special case of ^ = (parallel wave propagations relative to B). 
Dispersion relation {Bl) gives three exact solutions. 



t/' = 1 , (S9) 
C/' = n^ (SIO) 

For the modified Alfven wave speed = + W^/c^y/"^ greater (or less) 

than the thermal sound speed a, solution (-B8) corresponds to the fast (or slow) 
mode, while solution (-B9) corresponds to the slow (or fast) mode. Solution (SIO) 
represents the suprathermal mode. 

For the suprathermal mode (-B10) with uj'^ = k'^P, one gets p = 0, p = 0, 
= 0, Uy — 0, d ^ 0, P ^ 0, ^ 0, by = 0, and ~ 0. Physically, this 
is an acoustic wave in the suprathermal gas propagating along the background 
magnetic field B without disturbing the thermal gas. For a relativistically hot 
CRG, the value of b can be fairly close to the speed of light c. 

For Cd > o, the slow mode {B9) with lo^ — k^a^ gives p 0, p 7^ 0, 7^ 0, 
Uy = 5 = P = V;^ — by = 0, aud 6^ = 0. Physically, this is an 
acoustic wave in the thermal gas propagating along the background magnetic 

— * 

field B without disturbing the suprathermal gas. 

For Cjj > a, the fast mode {B8) with u"^ = k'^W'^ involves all nonzero per- 
turbation variables. By equation (-B3), p/po and S/S^ are out of phase, while by 
equation (-B5), p/po and b^/B are in phase as expected. 

For Cd < a, the fast mode {B9) with cj^ = /c^a^ has the same characteristics of 
an acoustic wave in the thermal gas propagating along the background magnetic 
field. One simply changes the name of this mode because of a faster wave speed. 

For Co < a, the slow mode {B8) with u'^ = k'^W'^ involves all nonzero per- 
turbation variables. By equation (-B3), p/po and S/6o are now in phase, while by 
equation (-B5), p/po and b^/B are now out of phase as expected. 

2. The case of ~ 7r/2 (nearly perpendicular wave propagations relative to B). 

Here, we would first correct typos in Parker's paper in the case of 9 being 
sufficiently close to 7r/2: Parker's solution (19) 

U^ = TT^ cos'^ (BU) 

M2 + 1 ^ ^ 

should be associated with the minus (not plus) sign of R for the slow mode, while 
his solution (20) 

+ 1 

^ cos' e J' \ (B12) 
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should be associated with the plus (not minus) sign of R for the fast mode, where 
R is exphcitly defined by his equation (18), namely 



R = ±[(M' + 1)' - 4M' cos' 9]^/' 



(BIS) 



(see his eqn. [15] or our eqn. [B7] with W^/c^ <C 1). Also, a more accurate 
version of his solution (21) for the suprathermal mode close to ^ = 7r/2 would be 



a 



(BU) 



as can be readily seen from his equation (15) by setting 9 = 7r/2. 

In the presence of displacement current effects, one can also derive three ap- 
proximate solutions from dispersion relation {B7) that may be written in the 
polynomial form of 



{n^ + 1) cos' 9 + M^ + 



(1 + a) sin' 9 



{l + a/n^){l + Wyc^) 



+ cos' 9 



cos' 9 + 



sin'^ 



n' + M'(n' + 1) 



(S15) 



U' - M'n' cos" 6* = , 



where = l^^/[a'(l-|-VF'/c')]. For 9 ~ 7r/2 and n' 1, the dispersion relation 
of the slow mode now becomes 



jj2 ^ 



M'cos'6' 



M2(l + l/n2) + sin' 9/(1 + + cos^ 9 ' 

the dispersion relation of the fast mode now becomes 

cos' 6'{[cos' 9 + sin' 6'/(l + M^'/c')]n' + M'(n' + 1)} 



t/' 



{v? + 1) cos2 ^ + M2 + (1 + a) sin ^/[(l + alv?){l + W^jc^)] 
and the dispersion relation of the suprathermal mode now becomes 



C/2 ^ M' + 



q; + 1 



+ 



am? cos' 9 



(l + a/n2)(l + W^Vc2) M2 + a + l ' 



(S16) 



(S17) 



(S18) 



For the slow mode (S16), and p/po are out of phase by using equation 

(-B5) as expected, while bjbo and pj are in phase by using equation (-B3). 

For the fast mode (-B17), h^jB and p/po are in phase by using equation (-B5), 
while 5/5o and p/po are out of phase by using equation (-B3). 

For the suprathermal mode (-B18), h^jB and p/po are in phase by using equa- 
tion (55), and bjbo and p/po are also in phase by using equation (53). 

3. The case of n' cos' 9 being sufiiciently larger than 1 with ^ 7^ 0. 

Corresponding to the suprathermal mode of Parker's solution (17), we now 



have 



ctsin 9 
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and one can readily show that 5/5o and p/po are in phase by equation (-B3), and 
hz/B and p/po are also in phase by equation (S5). 

Corresponding to the fast and slow modes of Parker's solution (16), we now 
have 



U 



2 [M^ + cos''e + si-ne/{l + W^/c^) + Rn] 



2 

r _ a[M^ + cos' e + sin' 6 /{I + W^c^) + Rp] tan' 9 \ 



+ 



1 



.n*cos4 6'. 

where, replacing definition (-B13) for R, R^ is defined by 

Rn = ±{[M' + cos' 9 + sin' 9/{l + W^'/c')]' - 4M' cos' ^}^/' (S21) 

with the plus and minus signs of Rd corresponding to the fast and slow modes, 
respectively. 

For the fast mode {B20) with the plus sign of R^ in definition (S21), one can 
show that d/do and p/po are out of phase by equation (S3), while b^/B and p/po 
are in phase by equation (-B5) as expected. 

Meanwhile, for the slow mode (-B20) with the minus sign of Rd in definition 
(521), one can show that S/S^ and p/po are in phase by equation {B3), while 
6^/5 and p/po are out of phase by equation {B5) as expected. 

Wc now summarize the results for the parameter regime of n ^ 1 (i.e., a) 
with comparable W and a. 

1. For the suprathermal mode, d/do, p/Po, and b^/B are all in phase. 

2. For the fast mode, 5/5o and b^^/B are out of phase with p/po and b^^/B being 
in phase. 

3. For the slow mode, 6/6^ and b^/B are out of phase with p/po and fe^/S being 
out of phase. 

APPENDIX C 

The mass conservation for the interstellar medium (ISM) is 

^^+V-{pv) = 0, {CI) 

where p is the ISM mass density and v is the bulk flow velocity of the ISM. After 
the operation of vertical integration / dz, one has 

^+V■{^^v) = 0, (C2) 
where p = / dzp is the surface mass density of the ISM. 
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Similarly, the mass conservation for the cosmic-ray gas (CRG) is 

^ + V-(5tI)=0, {C3) 

where 5 is the mass density and u is the bulk flow velocity of the CRG. After the 
operation of vertical integration / dz, one has 

— +v-(sii) = o, m 

where S = / dz5 is the surface mass density of the CRG. 

We invoke the polytropic approximations for both the ISM and the CRG, 

P = Clfi {C5) 

and 

P = C^S , (CQ) 

where C5 is the thermal sound speed of the ISM and Cc is the effective sound 
speed of the CRG, p is the vertically integrated (two-dimensional) gas pressure of 
the ISM, and P is the vertically integrated (two-dimensional) effective pressure of 
the CRG which is regarded as a "relativistically hot" gas. 

For a plasma of inflnite conductivity, the magnetic induction equation is 

^ = VxivxB), {CI) 

— * 

where the magnetic field B satisfies the divergence-free condition of 

V • 5 = . (C8) 

— * 

In the current context, bulk fiow velocities v±_ and u±_ perpendicular to B are 
constrained to be the same by the large-scale mean magnetic field. 

By ignoring effects of displacement current for nonrelativistic flows, we have 

VxB = —{Jism + Jcrg) , (C9) 
c 

where c is the speed of light, Jism and Jcrg are the electric current densities 
present in the ISM and in the CRG, respectively. 

The (two-dimensional) momentum equation for the ISM is 

/.ff + Vp- /d /^""''^ -/.V^ = 0, (CIO) 
Dt J c 

and the (two-dimensional) momentum equation for the CRG is 

Dt\mc^ J J c mc^ 

where m is the mean particle mass in the CRG, $ is the negative gravitational 
potential, i/j = mc^-|-rmP/[(r — 1)S] is the proper speciflc enthalpy including that 
of the particle rest mass. Given the polytropic approximation, i/j is a. constant. 

For the background rotational equilibrium in the radial direction, the above 
two equations can be combined to derive equation (2.1) (denoting e = 'Ei/j / [fhc^]) . 
Likewise, perturbations of these two equations can be combined to derive the 
radial component equation (2.6) by using equation (2.1) and the force-free ap- 
proximation. 



